One qubit almost completely reveals the dynamics of two 
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From the time dependence of states of one of them, the dynamics of two interacting qubits is determined 
to be one of two possibilities that differ only by a change of signs of parameters in the Hamiltonian. The 
only exception is a simple particular case where several parameters in the Hamiltonian are zero and one of the 
remaining nonzero parameters has no effect on the time dependence of states of the one qubit. The mean values 
that describe the initial state of the other qubit and of the correlations between the two qubits also are generally 
determined to within a change of signs by the time dependence of states of the one qubit, but with many more 
exceptions. An example demonstrates all the results. Feedback in the equations of motion that allows time 
dependence in a subsystem to determine the dynamics of the larger system can occur in both classical and 
quantum mechanics. The role of quantum mechanics here is just to identify qubits as the simplest objects to 
consider and specify the form that equations of motion for two interacting qubits can take. 



I. INTRODUCTION 

What can we learn about the dynamics of two interacting 
qubits by observing the time dependence of states of one of 
them? How much can the dynamics of the two be changed 
without changing the time dependence of states of the one? 
We will show that the dynamics of the two can be changed 
only by a change of signs of parameters in the Hamiltonian. 
There are only two simply related possibilities for the dynam- 
ics of the two qubits that give the same time dependence for 
the states of the one. This is set out in sections UlllI VI 

There is only one exception: in a simple particular case 
where several parameters in the Hamiltonian are zero, one of 
the remaining nonzero parameters can vary over the whole 
range of real numbers without changing the time dependence 
of states of the observed qubit. This is described in SectionlV] 

Determining the dynamics of the two qubits generally takes 
three time derivatives at time zero, expansion to third order in 
powers of time, of the mean values that describe the states of 
the observed qubit evolving in time, but it can take six time 
derivatives when some of the parameters in the Hamiltonian 
are zero. 

The mean values that describe the initial state of the unob- 
served qubit and of the correlations between the two qubits 
also are generally determined to within a change of signs by 
the time dependence of states of the observed qubit. This gen- 
erally takes four time derivatives. There are many exceptions. 
This is described in Section[yi] An example that demonstrates 
all the results is worked out in Section IVlIl 

Various implications and applications can be considered. 
One broad view of our results, framed only by the context of 
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open quantum dynamics of a system evolving together with its 
environment, is that time dependence in the system can pro- 
vide almost complete information about the environment and 
the interaction of the system with the environment, without 
measurements of the environment Our results provide 

an example that demonstrates the general statement. They 
have immediate application when the system is a qubit and 
its interaction with its environment can be modeled by inter- 
actions with qubits |3]. 

A more specific application can be seen in quantum infor- 
mation processing. When a physical device's performance of 
an operation is tested to verify that an interaction between two 
qubits is what it was designed to be, complete quantum pro- 
cess tomography fliSSiiElIiliiH is not needed. 
Our results suggest a simpler procedure. The dynamics of the 
two qubits can be determined almost completely by measure- 
ments of time dependence for one qubit with varying initial- 
ization of the state for that qubit but not for the two qubits. 

The time dependence of the states of the one qubit is an 
example of an open quantum dynamics described by maps of 
states of a subsystem caused by unitary Hamiltonian dynamics 
in a larger system. Since these maps generally do not depend 
on time as a semigroup, it is an open question how and to what 
extent a map at one time or at several times determines the 
map at other times. Here we give an answer for the dynamics 
of two qubits. The maps of the states of the one qubit for all 
times are determined by the Hamiltonian for the two qubits 
and the initial correlations between the two qubits 11141 lisll . 
We show how and to what extent these are determined by the 
maps in a neighborhood of the initial time. 

A classical analog, considered in Section lVIIII exhibits the 
logical structure of our method in a setting where it can be 
easily seen. Feedback in the equations of motion is what lets 
time dependence in a subsystem determine the dynamics of 
the larger system. This can happen in both classical and quan- 
tum mechanics. Here, quantum mechanics simply identifies 
qubits as the simplest objects to consider and specifies the 
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form that equations of motion for two interacting qubits can 
take. 



II. METHOD 

We consider two qubits, one described with Pauli matrices 
Si, S2, S3 and the other with PauU matrices Si, S2, S3. We 
assume the two qubits interact with each other but not with 
anything else. The time dependence is generated by a Hamil- 
tonian 

H = \oL3^i + ^/3feSfc + i7jfcSjSfe (2.1) 

with real parameters aj, Pk, Ijk- When an index j or k is 
repeated, there is to be a sum over the values 1, 2, 3 for that 
index. 

In almost all cases, the dynamics of the two qubits is almost 
completely determined by the time dependence of states of the 
S qubit alone. To see this, we look at the mean values 

(E„)(i) - (e^*«S„e-^*^) 

= u„j(t){T.j) + Vnkit){Ek) + w„jfc(SjSfe) (2.2) 

for n = 1 , 2, 3 which describe the states of the E qubit at time 
t. We consider variable initial states of the S qubit described 
by variable mean values (Ej). The u„j(t) are determined by 
the (S„)(<) for variable (Sj). We will see that the dynam- 
ics for the two qubits is almost completely determined by the 
Unj{t). Later we will see that generally the (Sfc) and (EjSfc) 
also are almost completely determined, so the initial state of 
the S qubit and of the correlations between the two qubits 
is almost completely determined by the time dependence of 
states of the S qubit, but nothing more can be learned about 
the dynamics from the (Sfc), (S^Sfc) and u„fc(t), Wnjk{t). 

The same dynamics may be described by different E^ , S^. 
There may be changes of the Ej, S^ that do not change the 
Hamiltonian, or that do change the way the Hamiltonian is 
expressed in terms of the Ej, Sfc but do not change the dy- 
namics. We will not be concerned with the differences made 
by these changes. 

The time dependence of states of the E qubit can be de- 
scribed 1I14I [isll by maps of the mean values (Ej). Of the 
(Ej), (Sfc), (EjSfc), the (Ej) describe the state of the E qubit, 
and the (Sfc) and (EjSfc) are considered to be parameters of 
the maps that describe how the dynamics of the two qubits 
drives the evolution of the E qubit ifTiifTsll . Different (Sfc) or 
(SjSfc) specify different maps. Each map applies to variable 
states of the E qubit described by variable (Ej). In almost 
all cases, the dynamics of the two qubits is almost completely 
determined by any one of these maps. 

III. MAGNITUDES 

We calculate Unj (t) as a series in powers of t by identifying 
the coefficients of Ej in the power series for e**^E„e~**^. In 
first order we get 

-i[T,i,H] = a2E3 -a3E2 +72fcE3Sfc -73fcE2Sfc (3.1) 



for n = 1. The coefficients of the Ej in this and the similar 
equations for — i[E2, i/] and — i[E3, i?] determine ai, a2, 

Ck3- 

In second order, the Ej terms in — « [ — « [Ei , H], H] are 

- [("2)^ + ("3)^ + 72fe72fc + 73fc73fe] El 

+ [aia2 + 7ife72fc]E2 + [aia3 + 7ifc73fc]E3. (3.2) 

These and the Ej terms in — i [ — i [E2, i/] , i?] and 
— i [ — i [E3, i/] , i?] determine 

7lfc7lfc, l2kl2k, J3k73k, llkl2k, 72fc73fc, 73fc7lfc (3-3) 

which are the dot products of the vectors 

71 = (711, 712, 713) 

72 — (721, 722, 723) 

73 = (731, 732, 733)- (3.4) 

Could anything more about the 7jfc be determined? The 
Hamiltonian is not changed if these three 7 vectors ( 13.4b and 
the vector {Pi, P2, Ps) are all changed by the same rotation 
when the Pauli matrices Si, S2, S3 are changed the same way. 
The dynamics is not changed for either qubit or for the system 
of two qubits by this change of the Si, S2, S3 used to de- 
scribe the S qubit. In particular, the mean values (E„)(i) are 
not changed. 

The sign of 71 • 72 x 73 is not determined, and it is not 
changed by rotations, so it is not changed when the Si, S2, S3 
are changed. For the same (E„)(i), there could be different 
possibilities for the dynamics of the two qubits corresponding 
to the two different possible signs of 71 • 72 x 73. We will con- 
sider these as two separate cases. We will find no equations 
that connect them. The sign of 71 • 72 x 73 can keep the same 
value in all the equations, either always plus or always minus. 
It is a free parameter It is not determined by the (E„)(t). 

For each sign of 71 -72 x 73, the 7jfc are determined to within 
rotations corresponding to changes of Si, S2, S3. We show 
in the Appendix that for each sign of 71 • 72 x 73, rotations of 
the Ej and Sfc can put the Hamiltonian in the form 

H = ^ajT,j + i/3fcSfc + i7fcEfcSfc (3.5) 

with real parameters aj, (3k, Ik- This change uses a rotation 
of the Ej and Unj{t) that can be found with knowledge of 
the dot products ( 13.31 ) that have been determined. For any 
dynamics of the two qubits, described by a Hamiltonian ( 12.11) . 
we can learn enough from the time dependence of states of 
the E qubit to change to a Hamiltonian of the form i3.5i and 
make the required change of of the Ej and Unj{t). From here 
on we will use this simpler Hamiltonian form to describe the 
time dependence. The changes made to put the Hamiltonian 
in the form ( 13. 5t do not change the value of 71 • 72 x 73, which 
gets called 717273. This also is shown in the Appendix. 

We will find that when 71 • 72 x 73 is zero, there still can be 
two possibilities for the dynamics of the two qubits that give 
the same time dependence for states of the E qubit. They dif- 
fer by a change of the signs of f32 and f]^. When the sign that 
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distinguishes the two possibilities is specified, the dynamics 
of the two qubits and the initial state of the S qubit and of the 
correlations between the two qubits are almost always com- 
pletely determined by the (I]„)(i). We have seen that ai, a^, 
aa are determined in first order and that (71)^, (72)^, (73)^ 
are determined in second order. 

In third order, in ~i\ , }i\ applied to Ei three times, the 
terms that involve things not already determined are 

7172/33S2 - 7371/52S3. (3.6) 

These and the similar terms in — i[ , applied to E2 and S3 
three times determine 7172/^3, 7273/?!, 737i/?2- We consider 
representative cases: 

[i) None of 71, 72, 73 are zero; 

(m) 7i 0, 72 7^ 0, 73 / 0; 
(«m) 71 = 0, 72 = 0, 73 7^ 0. 

In case (i), the magnitudes of (3^ are determined. In 

case (m), the magnitude of f3i is determined. In case (in), 
nothing new is determined. We will consider signs in Section 

In fourth order, in — i[ , H] applied to Ei four times, the Ej 
terms that involve things not already determined are 

(73)' [{f3if + m'] Si + (72)' [(/33)' + iPif] Si 

+71 72/3i/32i;2 +7173/31/33^3- (3.7) 

These and the similar terms in —i[ , H] applied to E2 and E3 
four times determine 

(73)'[(/3i)' + (/32)'], 7i72/3i/32, 

(72)'[(/33)' + (/3l)'], 7273/32/33, 

(7i)'[(/32)' + (/33)'], 737i/33/3i. 

In case (m), all the magnitudes of (3i, (32, /?3 are now deter- 
mined. In case (i) they were already determined. 

The rest of this section is for case {iii). There only (/3i)^ + 
{P2Y is determined. That is all that can be determined about 
/3i, /32. In case [Hi), the Hamiltonian contains /3i, /32 only 
in /3iSi + /32S2. It is not changed if (3i, (32 are changed as 
the components of a two-dimensional vector that is rotated 
when Si, S2 are changed the same way. The dynamics is not 
changed for either qubit or for the system of two qubits by 
this change of the Si, S2 used to describe the S qubit. In case 
(Hi), there are Si, S2 for which the dynamics is described by 
the Hamiltonian ( 13. 5t with (3i positive and (32 zero. Hence we 
assume that is the case. 

In sixth order, in — i[ , H] applied to Si six times, there 
is a term — (73)^(/3i)^(/33)^Si that, for case (iii), determines 
{(33)'^ when (/3i)^ is not zero. Case {iii) with /?i, /32 both zero 
is the one exception overall. It is described in SectionlV] 

IV. SIGNS 

All the magnitudes of the parameters a, (3, 7 are determined 
in all cases except the one described in Sec. |V] The signs of 



ai, a2, ct3 also are determined. We can change the signs of 
any two of Si, S2, S3. That allows us to choose Si, S2 S3 
so that neither 72 nor 73 is negative. Hence we assume that is 
the case. Then the sign of 71 is the sign of 71 • 72 x 73. The 
sign of (3i is determined in cases (i) and (ii) because 7273/31 
is determined. We are taking (3i to be positive in case (iii). 
Only the signs of 71, /32, (3^ are not determined. 

Can signs of 71, (32, (3^ be changed without changing the 
time dependence of states of the S qubit? In case (i), if the 
sign of any one of 71, /32, (3^ is changed, the signs of all three 
must be changed, because 7172/^3, 737i/32 and 7273/32/33 are 
determined. In case {ii), if the sign of (32 or (3j, is changed, 
the signs of both must be changed, because 7273/32/33 is de- 
termined. In case {iii), we are taking (32 to zero, so the sign of 
(3^ is the only one that can be changed. In all cases, the only 
change that can be made is the change of signs of all the 71, 
(32, (3z that are not zero. 

Changing the signs of 71, (32, (3^ does not change the time 
dependence of states of the S qubit. This change of signs 
relates two different Hamiltonians that describe different dy- 
namics for the two qubits but give the same time dependence 
for states of the S qubit. We will show this in two steps. In 
this section we show that the two Hamiltonians give the same 
Unj{t), because the change of signs makes no difference in 
any Sj terms in the power series for the e'*^S„e~**^. In 
Sec.|Vl]we will show that the two Hamiltonians give the same 
i;„fc(t)(Sfc) and Wnjkit){^j^k)- 

Let M be one of the Sj, S^, or SjS/j. Consider one of 
the times that AI occurs in a power series for a e'*^S„e~**''^. 
There are powers of the parameters a, (3, 7 multiplying M. 
Let p be the power of 71 plus the power of (32 plus the power 
of (3^. Now consider all the times that M occurs in the power 
series for the e'*^S„e~**''^. For each M, either p is even 
every time M occurs or p is odd every time M occurs. It is 
even for 

Sj, SjS2, SjS3, Si, (4.1) 

which we call blue operators, and odd for 

S,Si, S2, S3, (4.2) 

which we call red operators. To see how this happens, con- 
sider how p can change. The power series are generated by 
repeated application of [ , H]. Each [ , 71S1S1] brings in a 
power of 7i, each [ , /32S2] a power of (32, and each [ , /33S3] 
a power of /33, so p increases by 1 with each [ , G] where G 
is a term of H that is a parameter times a red operator, and p 
does not change when G is a parameter times a blue operator 
The commutator of two blue operators is a blue operator, the 
commutator of two red operators is a blue operator, and the 
commutator of a blue operator and a red operator is a red op- 
erator. The [ , G] that change p are the [ , G] that take blue 
operators to red operators and red operators to blue operators, 
and the [ , G] that do not change p are those that take blue 
to blue and red to red. A change of p between even and odd 
is a change of color When an M recurs, its p has changed 
between even and odd an even number of times. For each M, 
either p is even every time M occurs or p is odd every time M 
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occurs. Since p is zero for the E„ at the start, p must be even 
for the blue operators and odd for the red operators. 

In particular, p is even for the E^. The terms are not 
changed by the change of the signs of 71, (32, Ps- The Ej 
terms, the Unj{t), cannot distinguish the two possibilities for 
the dynamics of the two qubits. 



V. THE EXCEPTION 

The one exception is represented by case [iii) with f3i, (^2 
both zero. Then the Hamiltonian is 

= a,I],+73S]3S3+/33S3. (5.1) 

Here (3^ can vary over the whole range of real numbers with- 
out changing the time dependence of states of the E qubit; the 
/33S3 term commutes with the rest of H and with Ei, E2, E3. 



VI. TWO-QUBIT STATES 

The mean values (5fc) and (EjSfc) describe the initial state 
of the S qubit and of the correlations between the two qubits. 
They also are generally determined to within a change of signs 
by the time dependence of states of the E qubit. We have two 
possibilities for the dynamics of the two qubits. For each pos- 
sibility, the parameters a, (3, 7 are determined by the time de- 
pendence of states of the E qubit, so the Hamiltonian H and 
the Vnk{t), Wnjk{t) are determined. The mean values (E„}(t) 
that describe states of the E qubit in time provide linear equa- 
tions (|Z2] | for the (Sfe), (E^Sfc). The twelve (Sfc), (E^-Sfc) 
generally are determined by the four time derivatives in the 
power series for the (e**^E„e~**^) to fourth order for the 
three values of n. There are many exceptions. For example, 
in the case we have described lll4ll where 73 is the only one 
of the a, (3, 7 that is not zero, only (E1S3) and (E2S3} are 
determined, and when only 72 and 73 are nonzero, only (Si), 
(S1S2), (E1S3), (E2S3}, (E3S2) are determined; the other 
(Sfc) and (EjSfe) have no effect on the time dependence of 
states of the E qubit. 

The (2fe), (EjSfc) for one possibility for the dynamics of 
the two qubits are changed to the (Sfe), (EjSfc) for the other 
possibility by just changing the signs of the (EjSi), (S2), 
(S3). This follows from what we learned in Sec. |IV] Let M 
be one of the S^ or EjSfc. Every time {M) occurs in the 
power series for the equations ( 12. 2t it is multiplied by powers 
of the a, (3, 7 for which p is even if M is a blue operator, 
odd if M is a red operator Changing from one possibility for 
the dynamics to the other, changing the signs of 71, (32, (3^, 
changes the equations for the (S^), (EjSfe) by just changing 
the signs of the terms for which p is odd. That just changes the 
signs of the coefficients of the (M) for which M is a red oper- 
ator, the coefficients of the (EjSi), (S2), (S3). The equations 
for the (Sfc), (E^S^) for one possibility for the dynamics are 
the same as the equations for the other possibility for (S^), 
(EjSfc) with the signs of the (E^-Si), (S2), (S3) changed. 



Now we can see that the time dependence of states of the S 
qubit does not distinguish the two possibilities for the dynam- 
ics of the two qubits. The mean values (E„)(t) that describe 
the states of the E qubit in time provide Eqs. (12.2b . Some of 
the terms of the power series for these equations generally de- 
termine the (Sfc) and (E^-Sfc). When these (Sfe) and (EjSfc) 
are used in the remaining terms, the equations obtained are 
the same for the two possibilities for the dynamics of the two 
qubits. Again let M be one of the S^ or Ej-S^. If M is a blue 
operator, then both (M) and the powers of a, (3, 7 that multi- 
ply it in these equations are the same for the two possibilities. 
If M is a red operator, then both (M) and the powers of a, (3, 
7 that multiply it have the same magnitude and opposite signs 
for the two possibilities. The ability of Eqs. (12. 2t to distin- 
guish the two possibilities is not increased in the exceptional 
cases where some of the (S^) or (E^Sfc) are not determined. 

VII. EXAMPLE 

Here is a substantial example. Suppose ai, a2, and (3i, 
(32, (33 are zero. Then the dynamics generated by the Hamil- 
tonian (13.5b can be worked out very simply [15]. The three 
matrices EiSi, E2S2, E3S3 commute with each other (The 
different Ej anticommute and the different Sj anticommute, 
so the different EjSj commute). That allows us to easily com- 
pute 

(e"^Eie-"^) = (Sie-''^=^*^=^-^e~*^='*^^"^) 
— (El) cos72i cos73i 
+ (Si) sin72t sin73t 
-(E2S3) cos 72isin73t 
+ (E3S2) sin 72^008 73^ (7.1) 

using the algebra of Pauli matrices, and similarly 

(e**^E2e-^*^) = (E2)cos73icos7it 

+ (S2) sin73i sin7it 
-(E3S1) cos73tsin7it 
+ (EiS3)sin73tcos7ii, (7.2) 

{e'*"T,3e"''") = (E3)cos7itcos72t 

+ (S3) sin7ii sin72t 
-(E1S2) cos jitsmj2t 
+ (E2Si) sin7it cos72t. (7.3) 

These Eqs. (17. lb - ( 17. 3b give examples of the functions u„j (t), 
Vnk{t), Wnjk{t) in Eqs. (12. 2b . We consider the cases where 71, 
72, 73 are all nonzero. 

The way time dependence for the one qubit reveals the dy- 
namics of the two is particularly clear in these cases. The 
Unj (t) are even functions of t, so their power series have no 
first-order or third-order terms. This is not changed when the 
Ej and the Sfc are changed by rotations. Therefore, for any 
Ej and S^ that are used at the start, the results of the pro- 
cedures described in the first and sixth paragraphs of Section 
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|III]are that ai, ai, as are zero and that /3i, /?2, /^a are zero 
in the cases where 71, 72, 73 are all nonzero. In between, 
the procedure described in the second paragraph of Section 
Hm determines the dot products (13.3b . so the Hamiltonian can 
be put in the form ( 13.5b . The magnitudes (71)^, (72)^, (73)^ 
are determined. The signs of 71, 72, 73 are not determined. 
The Pauli matrices Si, S2, S3 can be chosen so that 72 and 73 
are positive. This leaves the sign of 71 undetermined. It is the 
sign of 71 -72 X 73. 

Thus, we see that time dependence for the one qubit reveals 
that the Hamiltonian can be put in the form ( 13. 5b and that ai, 
a2, 03 and /3i, /32, fij, are zero for the cases where 71, 72, 
73 are all nonzero. This gives Eqs. dT.lb - ( 17.3b . There are 
two different possibilities for the dynamics of the two qubits, 
corresponding to the two different signs of 71. We can also 
see that for each sign of 71, the Eqs.( l7.1b - ( 17.3b give the (Sfe) 
and the (E^Sfc) for j ^ k. The results for the two different 
signs of 71 differ in the signs of (S2S1}, (E3S1) and (S2), 
(S3), as described in SectionlVTl 

The three mean values (E^Sfc) are not determined. They 
are constants of the motion and they are not involved with the 
time dependence of the other mean values for the two qubits. 

In this example it is also particularly clear that there is a 
real difference between the two possibilities for the dynamics 
of the two qubits corresponding to the two different signs of 
71. There is no unitary transformation that changes one into 
the other. Changing the sign of 71 changes the eigenvalue 
spectrum of the Hamiltonian. The Hamiltonian (13.5b for these 
cases is a function 

if = i[-7l(S2S2)(E3S3)+72(S2S2)+73(S3S3)] (7.4) 

of the two matrices S2S2 and S3S3 which form a compete set 
of commuting operators. The four pairs of their eigenvalues 
label abasis of eigenvectors |1, 1), |1, —1), | — 1, 1), | — 1, —1} 
for the space of states of the two qubits. This shows that the 
eigenvalues of H are 

i[-7i +72 +73] 



^[71 
^[71 



■ 72 
72 ■ 



73] 
■73] 



-71 ~ 72 - 73j- 



(7.5) 



For positive 71, there is an eigenvalue of H that is less than 
any eigenvalue of H for negative 71. For negative 71, there is 
an eigenvalue of H that is greater than any eigenvalue of H 
for positive 71 . 



VIII. CLASSICAL ANALOG 

To see how much quantum mechanics is involved, we con- 
sider a classical analog. It exhibits the logical structure of our 
method in a setting where it can be easily seen. Let x and y' 
be real variables that have linear equations of motion 



dx 
'dt 



ax + J y 



dt 



with real parameters a, (3, 7', 5. We assume that neither 7' 
nor 5 is zero, and if 7' is negative we change the signs of 7', 5 
and y' to make 7' positive. Let 



y = 




7 = vyM- 



Then 



dx 



dy 
dt 



±jx + l3y 



(8.2) 



(8.3) 



with the ± the sign of S, and 



fx 
dt^ 
d^x 
'dfi 



a^x ± 7^x + 7(a + /3)y 

(a^ ±2a7^)x±7^/3a; 

+ {a^j + a(3j + (3'^-f±j^)y. (8.4) 



5x + /3y' 



(8.1) 



Looking at the time derivatives of x for variable x, we learn 
a from dx/dt, then 7^ and the ± sign from dx'^ /dt^, and f3 
from d^x/dt^. The parameters in the equations of motion for 
X and y are determined from the power series for the time 
dependence of x for variable initial values of x. The initial 
value of y is determined from the terms in the derivatives of 
X that do not depend on x. This is closely analogous to the 
calculations of Sec. Hill 

When neither 7' nor S is zero, there is a y for which the 
dynamics of x and y is described by the equations of motion 
( 18.3b . Then there is feedback from the equation of motion for 
y to the time dependence of x that depends on the initial value 
of X and allows all the parameters in the equations of motion 
for X and y to be determined from the time dependence of x 
alone. 

For qubits, we needed quantum mechanics just to identify 
qubits as the simplest objects to consider and to specify the 
form that equations of motion for two qubits can take. Every- 
thing we learned can be obtained from the equations of motion 
for the mean values (Sj), (Sfe), (SjSfe). There is no reason in 
principle that these could not be classical equations of motion 
for some system. They could even be Hamiltonian equations, 
with a Hamiltonian function of the form ( 13. 5b . for classical 
spin variables with suitably defined Poisson brackets Il6ll . 



IX. DISCUSSION 

We have shown that the Hamiltonian for two qubits can 
be determined almost completely from the time dependence 
of states of one of the qubits, but it requires rather detailed 
knowledge of that time dependence. In our calculations we 
generally need time derivatives up to third order, and in some 
cases up to sixth order. A practical application would be diffi- 
cult. 

The result seems to raise a question; it suggests that a step 
might remain to be taken to understand the mathematical re- 
sult physically. Why is the time dependence of states of the 
one qubit unchanged when the Hamiltonian for the two qubits 
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is changed in just those signs of three terms? Is there a rea- 
son, a physical explanation, perhaps a symmetry, that would 
let us predict this result without doing the calculations? Is 
there more to be said about that change of three signs? We 
have not found an answer. It may be that this result is just 
what happens mathematically in this particular situation, not 
a consequence of something more general, more physical, or 
more easily understood. 

APPENDIX: DIAGONALIZING THE GAMMAS 

Here is a simple proof that rotations can diagonalize the 
matrix of coefficients for the interaction terms and put 
the Hamiltonian in the form i3.5i . and that the rotation of the 
Sj and Unj (t) that can be found with knowledge of the dot 
products ( I3.3l l. The 7 vectors are central. Their dot products 
7m ■ 7„ are the elements of a real symmetric 3x3 matrix, so 
there is a real 3x3 rotation matrix R with elements Rjm such 
that 

Rjmim ■ Rkiiln = Rjm 7m ' In [R^^)nk = (A.l) 

for i ^ k. The three vectors Rjmlm are orthogonal. There 
is a rotation S that changes these three vectors and the Si, 
S2, S3 together, as described in Section Hill and takes each 
Rjmlm to a vector SRjmlm that is along the j axis. Then 

Ijk^jEk = SjT,- • S 

{Rjm^m)Rjn'~f7i ' ^ 
— {Rjm^m)SRjnln ' 
= iRjm^m)lj{SE)j 

= lj{Rjm^m)iSjlEi) (A.2) 



where is the only component of SRjmlm that may be 
nonzero, the component in the j direction, which may be pos- 
itive, negative, or zero, and the Sjk are the elements of the 
rotation matrix for S. 

The El, Yj2, S3 are rotated by R and the Si, S2, S3 are 
rotated by 5*. The rotation R depends only on the dot products 
Im'ln, which are determined by the time dependence of states 
of the Yj qubit. The 7 vectors are changed differently by the 
two rotations. The vector SRjmlm is just the vector Rjmlm 
rotated by S. The vector Rjmlm is not the vector 7,„ rotated 
by R; it is the linear combination with coefficients Rji, Rj2, 
Rj3 of the vectors 71, 72, 73. 

There is no change in the value of 71 -72 x 73, only a change 
of what it is called. It is 

717273 = (SRull) ■ {SR2mlm) X [S Rsnln) 
= {Rllil) ■ {R2mim) X (Ranin) 
= {R11R22R33 — R11R23R32 
+ R12R23R3I — R12R21R33 

+R13R21R32 — ^13-^22-^31)71 • 72 X 73 
= dct(i?)7i • 72 X 73 

= 71 • 72 X 73 (A.3) 

because det(i?) is 1 for a rotation. 
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